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The  Grant  commenced  on  1st  June,  1983,  initially  for  twelve  months, 
but  was  extended  for  a  further  year  to  31st  May,  1989.  The  research,  to  be 
carried  out  at  the  University  of  Southampton,  was  to  study  turbulent 
transport  in  confined  turbulent  flames.  Of  speeial  interest  was  the 
process  of  turbulent  transport  driven  by  a  pressure  gradient 
The  experimental  programme  was  .to  involve  the  simultaneous  use  of  Laser 
Doppler  Anemonetry  and  elastic  (Mle)  scattering  from  added  seed  particles 
to  evaluate  turbulent  mass  transport  fluxes  In  the  flame.  A  parellel 
theoretical  study  was  to  Improve  models  of  turbulent  transport  and  hence 
to  Identify  a  suitable  model  for  combustor  flow  field  calculations. 

However,  the  Grant  was  terminated,  at  the  request  of  the  University, 
with  effect  from  30th  November,  1984.  There  were  two  reasons  for  this. 

The  first  was  the  long  ten  Illness  and  eventual  resignation  of  the 
Research  Assistant  who  had'  been  appointed  to  carry  out  the  experimental 
research.  As  a  consequence  of  bis  aosence  essentially  no  progress  was  made 
In  the  experimental  programme.  The  second  reason  was  the  departure  of  the 
Principal  Investigator  to  take  up  an  appointment  at  the  University  of  Cambridge. 

This  Final  Report  Is  divided  into  three  parts.  The  first  deals  briefly 
with  the  abortive  experimental  programme,  while  the  second  Is  concerned  with 
theoretical  model  development  where  good  progress  was  made.  The  third  part 
of  the  report  briefly  discusses  the  question  of  models  for  combustor  flowfleld 
calculstlons. 

2.  Experiments 

The  experlmentst  rig  was  a  small  ducted  burner  with  a  rectangular  cross 
section  Incorporating  four  scblleren  quality  vitreous  fused  silica  windows 
and  supplied  with  premlxed  propane  and  air  at  velocities  up  to  20  m/s. 

The  flow  was  uniformly  seeded  with  T^O,  particles  which  permitted  simultaneous 
measurements  to  be  made  of  gas  density  and  axial  flow  velocity.  The  rig 
and  experimental  techniques  have  been  described  elsewhere^^ ^ . 

Although  the  Research  Assistant  learned  to  use  these  techniques  and 
successfully  reproduced  some  earlier  published  data,  no  new  data  of  publishable 
quality  was  obtained.  The  reason  for  this  failure  was  the  long  term  illness 
and  eventual  resignation  of  the  Research  Assistant. 


ExperlffleataX  data  obtained  by  these  nethods  is  described  in  earlier 
publications 

3.  Theoretical  Work 


The  objective  of  the  theoretical  profraaM  was  the  exploration  of 
theoretical  models  suitable  for  prediction  of  dtap  combustor  flowflelds. 

The  starting  point  of  this  work  was  the  model  of  prealxed  turbulent  combustion 
developed  by  the  Principal  Investigator  in  collaboration  with 
Professor  P.  A.  Libby  of  University  of  California,  San  Diego  and 
Professor  J.  B.  Moss  of  Crsnfield  Instiutte  of  Technology^^  ' 

Important  progress  wss  made  in  two  areas  during  tbs  period  of  tbs  grant 
as  described  below. 

3.1  Scales  of  Scalar  Field 


( i  > 

Central  features  of  the  model  of  premixed  turbulent  combustion  ’ 
which  forma  the  subject  of  this  research  are  (i)  the  description  of  the 
theraochemistry  in  terms  of  a  progress  variable  c  which  may  be  defined  as 
a  normalised  product  species  concentration  and  (ii)  consideration  of  a 
combustion  regime  in  which  the  turbulent  flame  consists  of  thin  moving 
interfaces  at  the  boundaries  between  packets  of  unburned  and  fully  burned  gas. 
The  combustion  reactions  then  occur  only  within  these  interfaces.  In  this 
combustion  regime  the  time  average  heat  release  rats  is  controlled  by  the 
length  or  time  scales  which  determine  the  interface  structure  and  its  frequency 
of  ocetirrence. 

fS\ 

The  instantaneous  equation  for  the  progress  variable  c  is  written 


3^(00  +  ^  (0U^C)  .  <p  D)_  —  ♦  w 
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where  w  is  the  chemical  source  term  and  subscript  a  represents  a  mean  value. 


For  statistically  stationary  conditions  the  mean  of  this  equation  at  high 
Reynolds  numbers  where  molecular  transport  nay  be  neglected  is 


where  w  is  the  time  average  of  the  source  term  w.  Our  objective  is  to 
derive  an  expression  for  w. 


(3) 


\ie  write 


Wj,(x)  \;(x) 


where  v(x)  '  Is  the  frequency  at  which  tlsMlet  interfaces  cross  location  x 


'r-</ 


w  d  t  > 


is  the  cheeieal  reaction  occurring  per  crossing,  wbeve  subscript  F  indicates 
integration  through  a  flaaelet  interface  and  the  triangular  brackets  denote 

a  tine  average. 

(*) 

The  analysis  nay  be  divided  into  three  parts.  In  the  first  part  it 
is  noted  that,  if  the  flaaelet  interfaces  are  thin,  c(x,  t)  at  a  fixed  point 
X  is  a  square  wave  tiae  series  between  c  >  0  and  c  >  1.  This  is  aodelled  as 
a  randoa  telegraph  signal,  leading  to  the  result  that  its  autocovariance 


Pll(*i  ■)  = 


<  C<£,  t)  C(i,  t  ♦  T)>  ■  cfl  -  c)  exp  [-  C* 

~  2c(l-c> 


The  crossing  frequency  v  is  deterained  as 


2c(l-c)  ,A 
V  a  '/T 


where  T  is  an  integral  tine  scale: 


r-  Pii  -  c* 

J  C  <l-c 


A  second  stage  of  the  analysis 


is  the  foraulation  of  a  balance 


equation  for  T.  This  is  done  by  writing  Equation  (1)  at  a  location  xi 
and  tine  ti  and  aultiplying  it  by  c  evaluated  at  location  and  tiae 
t, .  Averaging  then  leads  to  a  general  two-point  two-tiae  covariance  balance 
equation.  The  equation  is  specialised  to  one  point  and  two  times  when  it 
becomes  a  balance  equation  for  the  quantity  Pn  defined  in  Equation  (S). 
Closure  assuaptions  are  made  and  the  equation  is  integrated  in  t  to  give 
the  required  balance  equation  for  T(x).  The  balance  equation  shows  that, 
with  the  nodelling  assumptions  already  introduced,  T  is  a  constant. 

Main  features  of  the  above  analysis  have  been  confirmed  experimentally 

(8  •) 

as  reported  in  the  original  publication -and  since  verified  elsewhere 


1 


The  loet  part  of  the  analysis  Is  to  oodel  the  chealcal  tens  w 

(5)  ^ 

which  appears  In  (3).  This  Is  represented  as 

“r 

•f-  r - =  <*> 

:a  cos  6 

where  o  Is  the  density  of  unburned  alxture,  S,  Is  the  laminar  burning 
r  ^  — 

veloxlty  of  that  alxture,  u  Is  the  Favre  mean  flow  velocity,  and  9 

'X, 

is  the  average  angle  between  u  and  the  direction  noraal  to  a  flaaalet. 

(t) 

More  recently  ,  Equation  (8)  has  been  replaced  by  an  expression,  based  on 
an  enseable  of  strained  laalnar  flaaeleta,  which  penslts  detailed  chealcal 
kinetics  to  be  Included. 

An  analysis  of  scalar  length  scales  in  prealasd  turbulent  flaaes  has 

(9) 

also  been  reported 

The  advantages  of  the  models  described  here,  in  coaparlson  with  earlier 
aodels,  are  that  steps  In  the  analysis  can  be  tested  experiment ally  and 

<1 

that  detailed  chemical  kinetics  of  arbitrary  coaplexlty  can  be  Incorporated 
3.2  Generalised  Model  for  Turbulent  Transport  Fluxes 


The  aodel  considered  here 


Is  what  Is  known  as  a  second  order  closure 


model  in  which  a  separate  balance  equation  is  foraulated  and  solved  for 
each  second  order  covariance  tera  such  as  a  Reynolds  stress  component  su'^  u^ 
or  a  mass  flux  coaponent  o  u'j  c".  It  Is  necessary  to  go  to  this  level  of 
sophistication  because  the  pressure  gradient  transport  tens,  which  form  a 
major  focus  of  the  work,  appear  here. 

(1  It) 

Previous  work  '  concentrated  on  planar  oae-dlaensional  turbulent 
flaaes  separating  regions  of  unlfora  flow  In  unburned  and  fully  burned  gases. 
Thus  probleas  of  aodelling  nonunlfora  turbulent  flowflelds  In  these  unburned 
and  fully  burned  regions  were  avoided. 


The  work  reported  here 


proposes  a  general  three  dimensional  set 


of  second  order  closure  equations.  Unlike  the  equations  used  previously 
(1  <•) 

this  set  goes  over  properly  to  second  order  equations  In  the  constant 
density  low  Mach  niaber  flows  upstreia  and  downstreai  of  the  turbulent  flame. 

A  feature  of  second  order  closure  aodels  Is  that  third  order  covariance 
teimis  appear  and  must  be  represented  In  terns  of  second  order  quantities. 

Tn  thr  present  vsr!;  this  is  facilitated  by  the  Introduction  of  conditional 
velocities  and  turbulence  quantities  in  reactant  and  product  packets. 


-  0 


Thus  for  examp le  ^ ^  ^ 

P  u’’  u"  c"  /  p  =  h  (1-c)  (l-2c)(u.  -u  )  (a.  -  ) 

i  J  I  ^ 

■*•  u'  u*.  -  u/  u' 

1  JP  i  Jr  J 


(9) 


where  the  flame let  Interfaces  are  thin.  Here  subscripts  r  and  p  refer 

to  conditional  quantities  in  resctants  and  products,  respectively,  so  that 

for  example  u,  a  cu, . 

ip  1 

The  model  derived^^'^  from  these  considerations  is 


pu"  u  ”  c" 
p"  i<i-cf) 


p  u."  u/'u  " 

—2 _ J-_S. , 


(1  -  2h 


f-i -  J  ^  i  i 

p  c(l-c)  p  c(l-c)  p 


(10) 


'V  A# 
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(1-20 


pu."  c" 


Jcd-c)  pi(l-c)  ^  2f(l-c) 


li. 


♦  f, 


P  tl  **  11  **  p  u  **  c 
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P  u”  u”  P  tt"  c" 
« _ ^  _  J  b  i 

'V..  'jm  - 


p  c(l-  c) 


—  \  '\l 

pc(l  -  c) 


^  pu  "  u” 
K.  in 
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p  u  "  c” 

^ _ 

—  'Va  \ 

p  c(l-c) 


(11) 


where 


R 


')<  p  u".  u* 


c  k 

- Z 

e  p 


k  n 


J_ 

3*. 


p  u  "  u  ” 
r-i  J 


•a  ■  "  •  ‘li.’'  * 


and  e  ,  K.  and  K,  ,  are  constants;  it  and  e  ars  tha  tulrbulence 
s  ijo  Iji 

kinetic  energy  and  dissipation  function  respectively. 


Modelling  of  other  terns  in  the  second  order  balance  equations  has  been 
(10) 

discussed  elsewhere  The  resulting  set  of  modelled  balance  equations 

has  been  solved  for  a  simple  one  dimensional  flow^''*^  with  satisfactory 
rp«ii*tv . 
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4.  Reconunendatlons  and  Conclusions 


A  major  objective  of  the  research  prograane  was  to  propose  a 

theoretical  model  for  premlxed  turbuleat  combustloo  which  would  be  suitable 

for  use  In  dump  combustor  flowflelds.  Processes  of  transport and 

(4 ) 

turbulence  production  due  to  pressure  gradients  were  to  be  taken  Into 

(10) 

account.  This  objective  has  been  met.  The  transport  model  recommended 

Is  that  described  above  In  Section  3.2.  The  turbulent  chemistry  model  should 

(0) 

be  that  outlined  In  Section  3.1  as  modified  la  a  recent  publication  . 

The  modification  permits  the  Inclusion  of  detailed  chemical  kinetic  mechanisms 
of  arbitrary  complexity  through  the  Introduction  of  an  ensemble  of  strained 
iMlcar  flame  lets. 

The  experimental  objectives  of  the  programme  were  not  achieved. 

This  was  a  severe  disappointment  and  the  need  for  such  e:q>eriments  remains 
although  some  of  the  gaps  In  the  experimental  picture  have  since  been  filled, 
for  example  la  Refs.  (6,  7,  12,  13). 
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progress  variable 
constant,  see  Eq.  (11) 
diffusion  coefficient 
see  £q.  (10) 
see  Eq.  (11) 

7avre  aean  turbulence  kinetic  energjr 
constant,  see  Eq.  (10) 
constant,  see  Eq.  (10) 
autocovariance,  see  Eq.  (9) 

tiae 

laalnar  burning  velocity 
integral  tiae  scale,  see  Eq.  (7) 
i-th  velocity  coaponent 
chenical  source  tera  la  Eq.  (1) 
see  Eq.  (3) 
spatial  coordinate 


c 

J 

V 

p 


Subscripts 

l.J.k.B 

P 

r 

Other 

(") 

C) 

(  )" 


Favre  aean  viscous  dissipation 
aean  flaae  angle,  see  Eq.  (8) 
flsaelet  crossing  frequency 
density 
tiae  delay 

spatial  coordinates 

products 

reactants 

aean  value 
Favre  aean  value 


Favre  lluctuation 
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